## Replication code for Section 2 of the Supplmentary Information for:
##  Abrajano, Elmendorf, and Quinn.
##  "Labels vs. Pictures: Treatment Mode Effects in Experiments About
##   Discrimination." Political Analysis. 
##
## Marisa Abrajano, Christopher Elmendorf, and Kevin Quinn
## 9/24/2017
##



## load the long (stacked) data
load("./PhotosLong-ALL.Rda")
load("./WordsLong-ALL.Rda")

## subset to the first row of each matchup stack
mydata.p.long <- mydata.p.long[mydata.p.long$time == 1,]
mydata.w.long <- mydata.w.long[mydata.w.long$time == 1,]


mydata.p.long$choice.missing <- is.na(mydata.p.long$choice)
mydata.w.long$choice.missing <- is.na(mydata.w.long$choice)




####################################################################
## info for Table 2 of Supp Info
tab.p <- table(mydata.p.long$choice.missing)
tab.w <- table(mydata.w.long$choice.missing)

choice.prop.test <- prop.test(c(tab.p[2], tab.w[2]),
                              c(sum(tab.p), sum(tab.w)))
####################################################################





####################################################################
## info for Table 3 of Supp Info
tab.p.time <- table(mydata.p.long$choice.missing, mydata.p.long$matchup.time)
tab.w.time <- table(mydata.w.long$choice.missing, mydata.w.long$matchup.time)

choice.time.prop.test <- prop.test(c(tab.p.time[2,], tab.w.time[2,]),
                                   c(apply(tab.p.time, 2, sum),
                                     apply(tab.w.time, 2, sum)))
####################################################################





mydata.p.long$int.missing <- is.na(mydata.p.long$int.mot.high)
mydata.w.long$int.missing <- is.na(mydata.w.long$int.mot.high)

mydata.p.long$ext.missing <- is.na(mydata.p.long$ext.mot.high)
mydata.w.long$ext.missing <- is.na(mydata.w.long$ext.mot.high)

mydata.p.sub <- unique(mydata.p.long[,c("respondent.ID", "int.missing",
                                        "ext.missing")])
mydata.w.sub <- unique(mydata.w.long[,c("respondent.ID", "int.missing",
                                        "ext.missing")])



####################################################################
## info for Table 4 of Supp Info
tab.p.int <- table(mydata.p.sub$int.missing)
tab.w.int <- table(mydata.w.sub$int.missing)

int.mot.prop.test <- prop.test(c(tab.p.int[2], tab.w.int[2]),
                               c(sum(tab.p.int), sum(tab.w.int)))
####################################################################



####################################################################
## info for Table 5 of Supp Info
tab.p.ext <- table(mydata.p.sub$ext.missing)
tab.w.ext <- table(mydata.w.sub$ext.missing)

ext.mot.prop.test <- prop.test(c(tab.p.ext[2], tab.w.ext[2]),
                               c(sum(tab.p.ext), sum(tab.w.ext)))
####################################################################


